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TIME-DEPENDENT  LINEAR  DAE’S 
WITH  DISCONTINUOUS  INPUTS* 

BY 

Patrick  J.  Rabier  and  Werner  C.  Riieinboi.dt^ 

Abstract  Existence  snd  uniqueness  results  are  proved  for  initial  value  problems  associated 
with  linear,  time-varying,  differential-algebraic  equations.  The  right-hand  sides  are  chosen  In 
a  space  of  distributions  allowing  for  solutions  exhibiting  discontinuities  as  well  as  “impulses’* 

This  approach  also  provides  a  satisfactory  answer  to  the  problem  of  “inconsistent  initial 
conditions”  of  crucial  importance  for  the  physical  applications.  Furthermore,  our  theoretical 
results  yield  an  efficient  numerical  procedure  for  the  calculation  of  the  jump  and  impulse  of 
a  solution  at  a  point  of  discontinuity.  Numerical  examples  are  given. 

1.  Introduction. 

In  this  paper,  we  prove  existence  and  uniqueness  results  for  initial  value  problems  asso¬ 
ciated  with  differential-algebraic  equations  (DAE's)  in  R” 

(1.1)  Ai-hRi  =  6, 

where  A,  B  are  smooth  time-voTjnnj  linear  operators,  and  b  belongs  to  a  class  of  distribu¬ 
tions  with  values  in  R"  containing  the  functions  that  are  sir.  ,oth  in  ( — oo,  0)  and  |0,  oo)  and 
have  a  discontinuity  at  the  origin.  Such  discontinuities  on  the  right  side  occur  frequently 
in  physical  problems  modelled  by  DAE’s.  For  instance,  in  electrical  network  problems  a 
discontinuity  of  6  may  correspond  to  the  operation  of  a  switch  at  a  given  time. 

The  existence  and  uniqueness  theory  for  problems  of  the  form  (1.1)  with  smooth  b  (and 
consistent  initial  conditions)  is  now  well  understood,  see,  e.g.,  (CS7J,  (KuM92),  [RR93a]  and 
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the  references  giwn  there.  But  elementary  examples  show  that  th.‘  setting  of  distributions 
is  indispensable  for  handling  discontinuous  right-hand  sides.  For  example,  if  A  is  constant 
with  =  0.  i?  =  /,  and  5  =  bcH  where  bo  6  R"  and  H  is  the  Heaviside  function,  then 
the  solution  of  (1.1)  (in  this  case  unique;  no  initial  condition  needs  to  be  or  should  be 
prescribed)  is  i  =  boH  —  Abob  and  hence  involves  the  Dirac  delta  distribution. 

When  A  and  B  are  constant  and  6  exhibits  jumps,  Laplace  transform  methods  are 
available  to  find  the  solutions  of  (1.1),  but  the  problem  appears  to  remain  open  for  time 
dependent  coefficients  A,  B.  This  is  the  case  considered  in  this  paper.  Our  results  depend 
essentially  upon  our  recent  work  [RR93a|  on  a  reduction  procedure  that  transforms  the 
distribution  solutions  of  (1.1)  into  the  distribution  solutions  of  an  explicit  ODE. 

The  “consistency”  of  initial  conditions  represents  another  topic  of  considerable  theoret¬ 
ical  and  practical  importance  in  the  study  of  DAE’s.  As  is  well  known,  even  for  smooth 
6,  (1.1)  will  not  have  a  solution  starting  at  arbitrary  points  to  6  R".  Rather,  existence 
of  a  solution  in  the  classical  sense  requires  that  lo  satisfies  certain  constraints  called  the 
consistency  conditions.  On  the  other  hand,  suppose  that  the  physical  process  modelled 
by  (1.1)  starts  at  time  t  =  0,  and  that  for  <  <  0  the  state  variable  x(t)  has  evolved  in  a 
way  totally  unrelated  with  (1.1).  If  lim  x(l)  =  lo  exists  this  xo  represents  a  natural  data 
value  for  the  initial  condition  at  t  =  0.  But,  since  xn  has  no  reason  to  be  consistent  with 
(1.1)  at  f  =  0,  the  mathematical  theory  only  provides  that  (1.1)  has  no  solution  for  this 
choice  of  initial  condition,  which  is,  of  course,  a  physically  unacceptable  statement. 

It  turns  out  that  the  consistency  question  is  closely  related  to  the  problems  addressed 
here.  By  viewing  this  question  as  that  of  extending  a  known  state  x(f )  for  t  <  0  lo  a  solution 
of  (1.1)  for  (  >  0  via  a  solution  of  (1-1)  in  CP'fR))",  we  show  that  the  ambiguity  can  be 
resolved;  From  iq  =  ^lim  i(f)  we  find  that  a  unique,  computable  jump  to  a  consistent 
value  occurs  at  f  =  0.  Furthermore,  for  problems  with  index  u  >  2,  the  sudden  transition 
between  xo  and  the  consistent  initial  value  may’  also  create  a  (computable)  impulse;  that 
is,  a  linear  combination  of  S  and  its  derivatives.  Fbrther  evidence  that  our  solution  is 
the  correct  one  is  provided  by  showing  that  it  is  the  limit  of  the  classical  solutions  of 
the  problems  obtained  by  smoothing  out  the  right-hand  side  near  t  =  0.  These  results 


complement  in  various  ways  those  already  obtained  for  problems  with  constant  coefficients 
in  [VLKSl],  lCo82l,  or  [093]. 

Section  2  gives  a  brief  review  of  the  reduction  procedure  developed  in  |RR93aJ.  Initial 
value  problems  for  (1.1)  are  then  considered  in  Section  3  for  right-hand  sides  in  a  class  of 
distributions  which  is  a  close  relative  of  the  class  C^p  of  “impulsive-smooth"  distributions 
introduced  in  [HS83].  The  application  of  these  results  to  the  problem  of  inconsistent  initial 
conditions  is  discussed  in  Section  4,  and  some  straightforward  generalizations  are  presented 
in  Section  5.  Finally,  Section  6  presents  some  computational  results  that  illustrate  the 
resulting  algorithms. 

2.  Reduction  procedure  for  linear  DAE’s. 

Let  A,B  G  C°°(R;  £(R"))  and  set  Ao  =  A,  Bo  =  B.  The  reduction  procedure  developed 
in  (RR93a]  generates,  under  appropriate  conditions,  a  new  pair  (Aj+i ,  Bj+j )  from  the  pair 
of  coefficient  functions  {Aj,Bj),  j  >0.  More  precisely,  set  r_i  =  n  and  assume  that 
Aj,Bj  €  C^CR; £(R'''"' ))  for  some  integer  0  <  ry_i  <  n.  Moreover,  suppose  that 

(2.1)  rank  A,(<)  =  ry,  Vt  €  R, 
where  0  <  ry  <  ry_i  is  a  fixed  integer,  and  that 

(2.2)  rank  A,(t)®Bj(t)  =  r^.,,  Vt  £  R, 

where  Aj(t)®  By(f)  €  £(R'''-'  x  R’"'-' ,R'''-' )  is  defined  by  A,(f)®  Bj(f)(u,  v)  =  Aj(f)u  + 

Bj(t)v. 

Under  the  conditions  (2.1 )  and  (2.2),  it  is  shown  in  (RR93a]  that  the  following  mappings 
exist: 

(i)  P,  6  C'”(R;  £(R’'J-' ))  such  that  Py(f)  is  a  projection  onto  rge  Aj(t),  Vf  £  R. 

(ii)  Cj  £  C*{R;£(R''',R'^>->)  such  that  C^(t)  £  G^R*-'  ,ker(?,(f)Bj(f)),  Vf  £  R,  where 
Qj  =  l-P,. 
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(iii)  Dj  €  C*‘(R;£(R^<-‘,R'-'))  such  that  Dj(t)  e  GKrge  A,(t).R'>),  Vt  e  R. 

With  Cj  and  Dj  as  in  (ii)  and  (iii)  above,  we  define 

(2.3)  —  DfAjCj,  —  Dj{BjCj  +  AjC,), 

so  that  €  C“(R; £(R'> )). 

If  (2.1)  and  (2.2)  hold  for  every  index  J  >  0,  the  procedure  can  he  continued  indefinitely. 
At  the  stune  time,  since  the  sequence  r,  is  non-increasing,  there  is  a  smallest  integer  u  >  0 
such  that  By  (2.1),  we  then  have  A„(t)  €  G£(R''*'‘'),  V(  €  R,  and  further 

reductions  produce  pairs  {Aj,Bj),  equivalent  to  the  pair  [A^.B„)  in  a  sense  defined  in 
(RROSaj.  The  integer  :/  >  0  is  called  the  index  of  the  pair  (A,B),  and  it  can  be  shown 
that  1/  is  independent  of  the  specific  choices  of  Pj,  Cj  and  Dj,  0  <  j  <  i/  —  I  made  during 
the  process. 

Remark  2.1.  For  constant  A  and  B  it  can  be  shown  that  (A.B)  has  index  t/  for  some 
1/  >  0  if  and  only  if  the  matrix  pencil  AA  +  S  is  regular,  and  that  u  is  exactly  the  index  of 
the  matrix  pencil  XA  +  B.  □ 

From  now  on,  when  referring  to  the  pair  (A,B)  with  index  i/  >  0,  it  will  always  be 
implicitly  assumed  that  the  reduction  was  possible  up  to  and  including  step  i'  (and  hence 
beyond);  that  is,  for  the  time  being,  that  (2.1)  and  (2.2)  hold  for  0  <  ^  <  v  (and  hence  for 
J  >  t'  +  !)■ 

Suppose  now  that  the  pair  (A,B)  has  index  v,  and  consider  the  DAE  (1.1)  with  b  £ 
(■D'(R))".  The  condition  (2.2)  is  equivalent  to  the  invertibility  of  lAj(<)  @  Bj{t)]^  and 
hence  to  ker  Aj(t)nkerBj(t)  =  {0},  or,  equivalently,  to  the  invertibility  of  Aj{t)Aj{t)^  + 
Bj{t)Bj(t)^.  We  now  define  sequences  uo,-  --  ,Ui/-i  and  l>o,  .  ..  .  ^i-  of  distributions  as 
follows:  Set  bo  =  b  and,  generally,  if  6j,  0  <  ;  <  i'  is  known,  construct  Uj  by  multiplying 
the  distribution  bj  by  the  operator  Bj[AjAj  +  that  is. 


(2.4)  uj  =  BfiA.Aj  +  B,Bfr'b„ 


Moreover,  for  0  <  j  <  f  —  1,  define 

(2.5)  bj^i  =  Dj[b,  ~  Bu,  -  Ait,), 

and 

(2.6)  r„-,  =CaCi  -  C„-i  €  C‘*(R;£(R''-‘.R")) 

(2.7)  v,,-]  =  ua  +  Cqu,  +  C0C1U2  +  ■  ■  ■  +  Co  •  ■  •  C^_2U„_i  6  I'D  (R))". 

In  (RR93a|  it  is  shown  that  a  distribution  x  €  (■P'(R))"  solves  (1-1)  if  and  only  if  i  has 
the  form 


(2.8)  I  =  +  1;,-,, 

with  r^_i  and  u„_i  given  by  (2.6)  and  (2.7),  respectively,  and  €  (I>'(R))’'‘'"‘  is  a 
solution  of  the  ODE 

(2.9)  x^  +  A~'B^x^  =  A~'b^. 

Naturally,  the  equivzdence  between  (1.1)  and  (2.9)  via  (2.8)  is  true,  in  particular,  for 
classical  solutions;  that  is,  when  (say)  b  €  C°°(R;R").  Then  also  u,  and  b,,  defined  by 
(2.4)  and  (2.5),  respectively,  are  of  class  C°°,  and  so  is  in  (2.7).  In  this  CEise,  (2.8)  also 
transforms  initial  value  problems  for  (1.1)  into  initial  value  problems  for  (2.9).  In  fact,  i 
solves  (1.1)  under  the  initial  condition  i(<o)  =  lo  for  fixed  <0  6  R  if  and  only  if  xo  verifies 

(2.10)  lo  =  r„_i(<o)Ji,o  +  t'»-i(fo). 

for  some  x„g  €  R*^'-'  which,  of  course,  is  necessarily  unique  by  the  injectivity  of  r„_i(/u). 
Such  values  lo  are  called  consistent  with  the  DAE  (1.1)  at  to-  Evidently,  if  x  solves  (1.1) 
then  the  values  x(<)  €  R"  are  consistent  with  (1.1)  at  t,  Vt  €  R".  Moreover,  initial  value 
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problems  for  (1.1)  with  consistent  initial  values  at  the  given  point  to  have  a  unique  (C°“) 
solution  while  initial  value  problems  with  non-consistent  initial  values  have  no  (classical) 
solution. 

It  is  an  interesting  fact  that  the  condition  (2.1)  is  essentially  superfluous  if  A  and  B 
are  analytic  (see  (RR93a]),  partly  because  in  that  case  (2.1)  automatically  holds  with  Cj  = 
mu  rank  i4^(t),  except  perhaps  at  points  of  a  subset  Sj  consisting  only  of  isolated  points 
in  R.  For  <  £  Sj,  we  have  dim  rge  Aj{t)  <  rj,  but  it  turns  out  that  an  ‘extended  range”  of 
Aj(t),  denoted  by  ext  rge  Aj{t),  can  be  defined  with  the  properties  that  ext  rge  Aj(t)  D 
rge  Aj(f)  and  dim  ext  rge  Aj(t)  =  Tj,  V(  €  R  (and  hence  ext  rge  Aj{t)  =  rge  Aj(t),  Vt  6 
R  \5j).  This  allows  for  the  construction  of  parametrized  families  Pj,  Cj  and  Dj  as  before, 
except  that  “ext  rge  Aj(t)”  now  replaces  “rge  j4,(<)”  everywhere.  Thus,  assuming  only 
that  (2.2)  holds  for  all  indices  j,  we  can  still  construct  a  reduction  (d^+i ,  Bj+j )  of  ( Aj.Bj), 
and  the  index  i/  of  (A,B)  is  defined  as  before.  But  now.  we  have  only  A^{t)  6  GLIR”'-' ) 
for  t  €  R  \5„,  and  (1.1)  reduces  to  (2.9)  via  (2.8)  only  if  ^u(<)  €  GI(R’'--')  for  every 
t  €  R.  Thus,  in  this  case,  the  invertibility  of  A„{t)  for  all  t  is  no  longer  guaranteed  and 
must  be  assumed  independently.  As  a  result,  future  reference  to  pairs  (A,B)  of  index  v 
“with  invertible  A^(<)  for  every  t  €  R”  should  not  be  viewed  as  a  redundancy  but,  rather, 
as  a  reminder  that  condition  (2.1)  can  be  dropped  if  A  and  B  are  analytic  but  that  then 
invertibility  of  A^{t)  is  no  longer  guaranteed  to  hold  for  all  t. 

All  indicated  results  extend  verbatim  to  the  case  when  R  is  replaced  by  an  arbitrary 

O 

interval  J  where  distributions  in  J  are  now  understood  to  be  distributions  in  J.  If  J 
contains  one  of  its  endpoints,  initial  value  problems  for  (1.1)  with  a  consistent  initial  vahie 
at  that  endpoint  can  be  considered  in  the  classical  setting. 

3.  Initial  value  problems  with  discontinuous  right  sides. 

In  [HS83],  Hautus  and  Silverman  introduced  the  class  Cj^p  of  “impulsive-smooth”  distri¬ 
butions  in  [0,  oo).  We  first  need  a  strsughtforward  variant  of  this  concept  for  distributions 
in  R.  Throughout  the  remainder  of  this  presentation,  R*  denotes  R  \  {0}. 
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Definition  3.1.  The  distrihution  i  €  TTfR)  i.s  said  to  he  iinpuisive-smooth,  i  €  Ciiiip(R*) 
for  short,  if  there  are  functions  €  C®“(R)  such  that  x~,pH  —  tl'(  \  —  H )  is  a  distribution 
ivith  support  {0}  where  H  denotes  the  Heaviside  function. 

If  X  e  Cimp(R*)  and  €  C°“(R)  are  such  that  i  -  <fi,H  -  0,(1  -  hf)  is  a 

distribution  with  support  {0},i  =  1,2,  then  (ipi —(,P2)-ff +  (</’!  ~t/’2)(l  —  ff)  is  a  distrihution 
with  support  {0}  and  hence  must  be  a  linear  combination  of  the  Dirac  6  and  its  derivatives. 
But,  since  it  is  also  a  function,  it  must  be  0;  that  \s,ipiH  +  tpi{\  —  H)  —  tpiH  +  0j(l  —  H). 
This  shows  that  i  —  ip,H  —  0,(1  —  H)  is  identical  for  i  =  1  mid  i  =  2,  and  hence  can  be 
called  the  impulsive  part  of  i,  denoted  by  limp- 

Therefore,  for  given  i  c  Cimp(R*)  the  difference  i  —  ijmp  has  the  form  tpH  +  0(1  —  hf) 
with  <0,0  €  C“(R).  Of  course,  <0  and  0  are  not  uniquely  determined  by  this  condition, 
but  0|(.„  ,)  and  <0|ju  are.  Thus,  there  is  no  ambiguity  in  setting 

~  ’^l(-oo.O)’  “  V’ljo  o„)- 

With  this  definition,  i_  €  C°^((-oo,0])  and  i+  €  C'~((0,oo)),  and  extending  x_  by  0  for 
t  >  0  and  1+  by  0  for  <  <  0,  we  may  write 

(3.1)  X  =  X_  +  X+  +  Xin,p, 

where  each  of  the  three  terms  on  the  right  side  is  uniquely  determined  by  x.  Conversely, 
given  x_  €  C°°((-oc,0])  and  i+  €  C'”((0,<xi))  and  a  distribution  ii^p  with  support  {0}, 
equation  (3.1)  defines  an  element  i  of  CimpfR*)- 

Remark  3.1.  Despite  the  terminology  “impulsive-smooth”,  it  should  be  kept  in  mind 
that  for  I  G  Cjn,p(R*),  x  —  Xjmp  is  not  a  smooth  function  in  R  since  it  may  have  a  jump  at 
0.  But  its  restrictions  x_  and  x+  to  (— oo,0)  and  (0,oo),  respectively,  extend  as  smooth 
functions  in  (— oo,0]  and  [0,  oo),  respectively.  □ 

Three  trivial  but  essential  properties  of  impulsive-smooth  distributions  are  the  following: 
(i)  Every  x  €  Cimp(R*)  may  be  a.ssigned  a  value  at  every  point  1^0,  namely  x(<)  =  x_(<) 
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if  <  <  0  and  i(t)  =  x+(0  if  f  >  0.  (ii)  The  derivative  and  the  primitives  (in  the  sense  of 
distributions)  of  i  €  CmpCR*)  are  themselves  in  Cimp(R').  (iii)  CimpCR')  is  both  a  vector 
space  over  R  and  a  C°°(R)-module.  In  fact,  if  i  €  CmpfR*)  and  (3,1)  is  used,  then  we 
have 


k  k->  /  J.  \ 

(3.2)  oi  =  ai_+ai+  +  53(^(-l)M^'^')  o<^*(0)A.+,)i'-\  Vo  6  C*(R), 

1=0  >=o  \  1  / 

whenever  limp  =  H!=o  A,  €  R,  0  <  «  <  Ir.  The  above  properties,  including  (3.1), 

have  an  immediate  generalization  to  elements  of  Ci^pfR*)  s  (Cimp(R*))".  In  particular, 
for  X  €  C;;„p(R*),  and  M  €  C^IR;  AR'.R")),  we  have  Mr  €  C;;;,p(R‘)  and 

t  t-i  ,.  , 

(3.3)  A/x  =  A/i_+Mi+  +  52(II(-lKr  j  A/'^>(0)A.+j)  A'’>, 

1=50  JSO  \  J  / 

whenever 

k 

(3.4)  iin,p=  A.  €R",  0<i<fc. 

1=0 


Deflnitioii  3.2.  Let  x  €  Cimp(R*),  so  that  there  exists  a  unique  decomposition  (3.4).  The 
impulse  order  of  x,  denoted  by  iord(x),  is  defined  as  follows: 

(i)  If  Xi  =  0,  0  <  i  <  k  and  i+(0)  =  i-(0),  and  hence  i  e  C”(R;R")),  then  set 


iord(x)  =  — m  —  2, 


where  0  <  m  <  oo  is  the  largest  integer  such  that  x  €  C’"(R;R"). 

(ii)  If  X,  =  0,  0  <  i  <  k  and  x+(0)  ^  X-(0)  (and  hence  x  has  a  discontinuity  at  the  origin), 
then  set 

iord(x)  =  -1. 

(iii)  If  X,  ^  0  for  some  0  <  i  <  k  then  set 

iord(i)  =  max{i  ;  0  <  i  <  h,  Aj  ^  0}. 
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Remark  3.2:  Let  M  e  C°°(R;£(R",R'")).  By  (3.3)  we  have  iord(A/i)  <  iord(x)  and 
equality  holds  if  n  =  m  and  Af(0)  is  invertible.  □ 

The  following  lemma  j>rovides  some  more  precise  preliminary  results  about  the  primi 
tives  in  the  sense  of  distributions  of  the  elements  of  i’iJ„p(R’)' 

Lemma  3.1.  (i)  Let  f  €  Cimp(®*)  have  impulse  order  k  €  Z  U  {  -oc}  and  let  y  £  (I>'(R))" 
be  such  that  y  —  f.  Then,  y  6  Cimp(R*)  and  !/  has  impulse  order  1-1 

(ii)  Let  f  e  f^jmp(R’)  and  let  Xg  €  R"  and  to  €  R*  be  given.  Then,  there  is  a  unique 
y  £  Cj^p(R*)  such  that  y  —  f  and  exactly  one  of  the  following  conditions  holds 

(3.5)  («)  y(<o)  =  io.  ((>+)  y+(0)  =  lo,  (i>-)  y-(0)  =  Jo. 

(Hi)  Let  the  sequence  /'  6  (!)i’J„p(R'),  f  >  1,  and  xq  €  R"  and  tg  6  R"  be  given.  Suppose 
that  there  are  an  open  interval  he  about  tg  and  some  f  €  Ci'Jnp(R')  such  that 

(36)  /V=K’ 

as  distributions  in  /(„  (i.c.  as  functions  ifO  ^  /,,)  and  that 

(3.7)  lim  /'  =  /  in  (P'(R))"- 

(—OO 

Let  i/  €  Cj'J„p(R*)  and  y  £  Cj';;p(R*)  be  such  that  y’  =  /',  y'(<o)  =  and  y  =  /, 
y(lo)  =  xg  (see  (ii)  above).  Then,  we  have 

(3.8)  limy'  =  y  in  (P'{R))". 

(A  similar  result  holds  if  to  =  0  and  y^,  y  are  char  cterized  by  y^  =  y4(0)  =  xg  and 

y  =  f,  y±(0)  =  Xg.) 

Proof:  (i)  Any  two  primitives  (in  the  sense  of  distributions)  of  an  element  of  (P'(R)J’’  differ 
from  a  constant  vector  of  R",  and  addition  of  a  constant  vector  dors  not  affect  membership 
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to  CtJnplR*)  nor  the  impulse  order.  Thus,  it  suffices  to  show  that  /  has  one  primitive  in 
C|'J„p(R*)  with  impulse  order  k  -  \ 

Write  /  =  /4.  +  /_  +  /,mp  with  /,„p  =  ^  *"•  choose 

<0  6  R.  Since  the  function  /+  +  /-  is  locally  integrable,  set  y  =  //^(/+  +  /-  )(s)ds.  so  that 

y  6  C“°((-oo,0J;R")  n  C*’(|0,oo);R")nC''(R;R")nC,"„p(R*)  and  y  is  a  primitive  of 
/+  +  /_  in  the  sense  of  distributions.  Furthermore,  it  is  obvious  that  y  €  C'"'*'' (R;  R" ) 
whenever  /++/_€  C'"(R;R"),  0  <  m  <  no.  Thus,  iord(y)  =  iord(/>.  +  /_ )  -  1  <  -2.  In 
pturticultir,  y  is  a  primitive  of  /  with  impulse  order  t  —  1  when  k  <  0  since  /  =  /+  +  /_  in 
this  case. 

Suppose  now  that  fc  >  0  and  set  yimp  =  X)!=i  that  /loH  +  y.^p  is  a  primitive 

"f /imp-  Evidently,  poW +  yimp  €  C|';„p(R')  and  iord(poff +  yimp)  =  1  >  -1  since  t  >  0 

and  fn  ^  0.  Thus,  y  =  y  +  noH  -t-  yimp  €  fimplR*)  is  a  primitive  of  f  ami  iord(y)  =  1:  —  1 
since  iord(t/)  <  -2  <  -1  <  iord(po/f  +  yimp). 

(ii)  The  primitive  of  /  obtained  in  (i)  verifies  y.  =  y_  and  y+  =  y+  +  po.  Since 

y  is  continuous  and  y(to)  =  0,  this  yields  y-(fo)  =  0  and  y+ffo)  =  Mo-  Every  other 

primitive  of  /  (still  denoted  by  y  for  simplicity  of  notation)  is  uniquely  characteri;:ed  by 
a  vector  Ag  e  R"  and  verifies  y-(to)  =  -^oiy-s(to)  =  Mo  +  -^o-  As  a  result,  Ag  =  ig  (resp. 
Ao  =  lo  -  Mo)  is  the  only  possible  choice  yielding  y_(fo)  =  ig  (resp.  y+(to)  =  ig).  Letting 
<0  =  0,  we  obtain  existence  and  uniqueness  of  y  €  Cjmpt^*)  such  that  y  =  J  and  either 
(3.5)(b+)  or  (5.5)(b_)  holds.  Next,  letting  tg  ^  0  and  observing  that  y(to)  =  y_(fg)  if 
<g  <  0  and  y(tg)  =  y+{to)  if  <o  >  0,  we  obtain  existence  and  uniqueness  of  y  €  r,’J„p(R’) 
such  that  y  =  /  and  y(/g)  =  ig. 

(iii)  To  begin  with,  let  us  briefly  recall  how  primitives  of  distributions  are  defined:  Let 
S  e  T>(R)  be  such  that  J^9  =  1.  For  y?  €  (I>(R))",  there  is  a  unique  £  (P(R))"  such 
that  V'  =  ip  -  S  y>  and  the  correspondence  yi  >-»  r/>  is  continuous  for  the  usual  topology  of 
(©(R))".  Note  also  that  supp C  suppi^U  supp  6.  Given  T  =  (Tj ,-- - 

the  formula 

(3.9)  (S,y»)  =  -(r,V')+c- 
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T 


with  c  €  R"  and  the  dot  denoting  the  usual  inner  product  of  R" ,  defines  S  as  a  distribution 
with  values  in  R"  and  shows  that  S  —  T,  and  all  the  primitives  of  T  are  of  the  form  (3.9) 
for  some  c  6  R" . 

In  general,  the  formula  (3  9)  does  not  permit  us  to  assign  a  value  S{t)  to  S  for  any 
t  €  R.  But  suppose  that  there  are  to  6  R  and  an  open  interval  /(„  about  to  such  that 
is  (say)  a  C”  function,  whence  T ■  v?  = /j  I"  for  all  v?  fe  W’e  may 

choose  6  such  that  supp  8  C  and  then,  for  t  £  we  may  define  5o(l)  =  /,„  7’(s)(fs; 

that  is,  S  =  (Soi,  ■  •  ,  Son)  with  So,{t)  =  T,(s)<fs.  In  (3  9)  let  c  =  (C], .  .  ,r„  )  with 

(3.10)  c,-J  So,(t)«(t)dl,  1  <  1  <  n. 

This  makes  sense  since  supp  6  C  and  So,(f)  is  defined  for  f  6  Let  tp  € 

whence  xl<  6  {t>{  '’(„))"  since  supp  8,  supp  <p  C  /i„-  As  T  =  So  in  (I>'(/|„))".  relation  (3.9) 

reads 


{S,<p)  =  -(So,V>) 

=  (So,t/0  +  c  ■  jf  s?  =  (So.yt  -  8  +  c  ■ 

=  {So,‘fi)  -  {SoJ-J  ip)  +  c 


But 


(So,8j^p)  =  J^ 

^  ^ I  I 

by  definition  of  c  in  (3.10).  Thus,  (S,ip)  =  {So-,p),  for  all  p  £  (I>(/i„))",  i.e.,  S|,_^  =  So. 
Decause  So  is  a  function  and  vanishes  at  <o.  it  follows  that  S  in  (3.9)  may  be  referred  to  as 
the  primitive  of  T  vanishing  at  to  when  r  is  chosen  as  in  (3.10)  (and  8  verifies  supp  8  C  /,J. 
The  iadependei-c  “  this  definition  from  the  choice  of  9  is  easily  seen:  If  S,  5  corresponds 
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to  two  such  choices,  we  have  5  =  5  +  c,c  6  R"  since  both  S  and  5  are  primitives  of  T, 
and  f  =  0  from  S|,^  =  So  =  S|,^  ,  As  a  result,  given  ig  €  K",  S  +  lo  may  be  referred  to 
as  the  primitive  of  T  verifying  S(to)  =  To- 

Now,  with  T  as  above,  let  e  (X>'(R))“,  f  >  1,  be  a  sequence  such  that  Tf  =  Ti, 
This  assumption  ensures  that  T*,  f  >  1,  as  well  as  T  define  the  same  vector  r  in  (.3.10). 
With  this  choice  of  c,  the  distribution  S*  6  (1?'(R))”  obtained  by  replacing  T  by  in  (3.9) 
is  the  primitive  of  T*  vanishing  at  to.  and,  under  the  assumption  lim  =  T  in  ('D'(R)}’’. 

i-^ao 

it  is  then  obvious  that  for  €  (P(R))"  we  have  lim  i.e.  lim  S^  =  S  in 

i—»oc  t^OO 

{V{R))" .  In  turn,  this  implies  that  lim  S*  +  lo  =  S  +  lo  td  6  R"  ■ 

/— •<» 

It  should  be  clear  that  part  (iii)  of  the  lemma  follows  from  the  above  considerations 
with  T  =  f  and  T*  =  /^,  so  that  S  +  Xa  =  tj  and  +  lo  =  y^-  The  proof  that  a  result 
similar  to  (3.S)  holds  when  <o  =  0  and  y*,  y  are  characterized  by  y^  =  f’.  y±(0)  =  lo 
and  y  =  /,  y±(0)  =  lo,  easily  follows  from  the  above  considerations  with  T  =  /++/_, 
=  /+  +  /i,  and  the  remark  that  =  /itnp.  Ihr  all  f  >  1,  since  and  /  coincide  as 
distributions  in  some  open  interval  about  the  origin  by  hypothesis.  Details  are  left  to  the 
reader.  □ 

Remark  3.3:  Because  of  condition  (3.6),  Lemma  3.1  (iii)  gives  an  unusual  result  about 
continuous  dependence  for  initial  value  problems.  The  incorrectness  of  this  result  under 
condition  (3.7)  alone  can  be  seen  even  in  the  case  when  n  =  1  tmd  f,f'  €  C*(R).  In  fact, 
by  the  theory  of  Fourier  series,  the  sequences  t^costt  and  (°  sin  (t  tend  to  0  in  I>'(R)  for 
every  a  €  R.  In  particular,  if  /  =  0,  /*(<)  =  (sin(i,  we  have  lim  /^  =  /  in  X>'(R)  as  in 

f-»oo 

(3.7).  Choosing  <o  =  0,  we  find  y  =  0,  y*(<)  =  1  —  cos  ft  in  Lemma  3.1  (iii)  But  then, 
lim  y^  =  1  /  y  in  I>'(R)  and  (3.8)  fails  to  hold.  □ 

t—oo 

Lemma  3.1  has  a  direct  application  to  initial  value  problems  for  the  ODE 
(3.11)  i  +  Mx  =  f, 

considered  in  the  following  theorem: 
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Theorem  3.1.  Let  A/  t  C^(R;£(R"))  and  let  /  €  C,'J„p(R*)  have  impulse  order  k  6 
Z  U  {-oo).  Then 

(i)  The  solutions  i  €  (I>'(R))''  of  the  ODE  (3.11)  belong  to  C/JnpCR*)  and  have  impulse 
order  it  —  1 . 

(ii)  For  given  Xo  6  R"  and  tg  €  R*.  the  ODE  (3  11)  together  with  one  of  the  initiaj 
conditions 


(3.12)  (a)  i(to)  =  io,  (fc+)  i+(0)  =  io,  (6_)  i_(0)  =  jo, 


has  a  unique  solution  x  €  C|'J„p(R*),  but  in  the  cases  (3.1 2)(b+)  and  (3.12)(b-  )  the  solutions 
corresponding  to  i+(0)  =  xo  and  i_(0)  =  to  need  not  be  the  same. 

(Hi)  Let  the  sequence  €  C/J„p(R*),  ^  >  1,  and  Xg  €  R"  and  to  €  R’  be  given.  Suppose 
that  there  is  an  open  interval  about  to  such  that 


(3.13) 


K 


V(>1. 


as  distributions  in  I,,,  (i.e.,  as  functions  ifO  ^  Ug)  and 


(3.14)  lim/'  =  f  in  (I>'(R))". 

^-•OO 

Let  X  (resp.  x‘)  €  C/J„p(R*)  denote  the  unique  solution  of  the  ODE  (3.11)  (resp.  (3.11) 
with  f  replaced  by  f*)  verifying  i(to)  =  *o  (resp.  i*(to)  =  xo)  whose  existence  is  ensured 
by  part  (ii)  of  the  theorem.  Then 


(3.15)  limx'  =  i  in  (P'(R))". 

1  —  00 

(A  similar  result  holds  if  to  =  0  and  the  initiaj  condition  for  x',  x  is  chosen  as  x±(0)  =  xg 
and  x±(0)  =  Xg,  respectively.) 


Proof:  (i)  Fix  to  €  R  and  denote  by  U  ^  C””(R;£(R”))  the  solution  of  the  initial  value 
problem  U  +  MU  =  0.  U{tg)  =  I.  It  is  well  known  that  t’(t)  C  G£(R")  for  all  t  €  R". 


Then,  z  €  solves  (3.11)  if  and  only  if  y  =  U~'x  6  (lytR))"  solve.s  the  equation 

y  =  U-'  f.  Since  /  €  Clp(R*),  we  have  V'f  €  Clp(R*)  and  iordfF-'/)  =  iord(/)  =  k 
from  Remark  3.2.  Next,  by  Lemma  3.1  (i)  the  solutions  of  y  =  U~'  f  are  in  C’.'J^plR*)  and 
iord(i)  =  iord(y)  =  t  —  1  by  another  application  of  Remark  3.2. 

(ii)  Since  [/(to)  =  /,  we  have  i(fo)  =  lo  (resp.  4  =  0  and  ii.(0)  =  id)  if  and  only  if 
y(fo)  =  lo  (resp.  to  =  0  and  y±(0)  =  Xo).  Existence  and  uniqueness  of  i  thus  follows  from 
existence  and  uniqueness  of  y  ensured  by  Lemma  3.1  (ii). 

(iii)  Prom  conditions  (3.13)  and  (3.14)  and  using  the  continuity  of  multiplication  of 

distributions  by  C*  matrix-valued  functions,  we  infer  that  U~'ff  =  and 

■'t© 

lim  U~'f*  =  U~'-f  in  (I5'(R))".  Denoting  by  y'  the  solution  of  y'  =  U~'  f*.  y'(0)  =  lo, 

f—OO 

we  find  that  lim  y*  =  y  in  (i>'(R))"  by  Lemma  3.1  (iii).  Thus,  lim  x'  =  i  in  (I^IR))" 

t—oo  t— oo 

since  x*  =  l/y*,  x  =  Uy  and  multiplication  by  V  is  continuous.  □ 

Remark  3.4:  Let  /i„,p  =  ^  0  From  Theorem  3.1,  we  have  ij„p  = 

EfJo  solution  i  €  Ci^p(R*)  of  i  +  Mi  =  /.  Comparing  impulsive  parts 

and  using  (3.3)  -  (3.4)  yields 

■  .^t-i  =  ^‘k^ 

£'(-1)'(-'|')m<»(0)A.+,  =  P..  1<.<1:-1. 

fc-l 

x+(0)  -  x-(0)  +  53(-1)>M<»(0)A,  =  pD. 

ISO 

By  inverting  these  formulas,  we  find  Aq,  •  •  •  ,  At_,  (depending  only  upon  Pi ,  •  •  . M*)  as 
well  as  i+(0)  -  i-(0).  Thus,  both  li^p  and  i+(0)  -  i-(0)  are  calculable  and  depend 
solely  upon  /imp.  □ 

We  now  focus  on  initial  value  problems  for  the  DAE 
(3.16)  Ax  +  Bi  =  b. 

w'here  A,B  €  C'*’(R;  £(R"))  and  6  e  C/J^pfR*)-  Under  the  assumption  that  the  pair  (A,B) 
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has  index  x/  >  0  in  R  and  that  A„(t)  is  invertible  for  every  <  €  R  (see  Section  2),  the  DAE’s 

(3.17_)  A{t)i.  +  B{t)x^  =  h^{t)  in(-oc,0), 

and 

(3.17+)  >l(t)i+ +  B(<)i+ =  6+(t)  in(0,oo), 

have  coefficients  and  right-hand  sides  of  class  C“°  in  (— oo.O]  and  [0,  oo),  respectively.  As  a 
consequence,  it  makes  sense  to  speak  of  values  lo  €  R"  which  are  consistent  with  (3.17_) 
(resp.  (3.17+))  at  a  point  to  <  0  (resp.  to  >  0)  in  the  sense  of  Section  2. 

Theorem  3.2.  Let  the  pair  (A,B),  A,B  €  C°°(R;£(R"))  have  index  v  >  0,  and  with 
the  notation  of  Section  2,  suppose  that  A„(t)  is  invertible  for  every  f  6  R.  ff  i)  6  C|^p(R*) 
has  impulse  order  h  €  ZU  {-oc},  then,  the  solutions  x  6  (I>'(R))"  of  the  DAE  (3.16)  are 
in  C|'[^p(R*)  and  have  impulse  order  at  most  k  +  i/  —  1.  Moreover, 

(i)  if  to  <  0  (resp.  to  >  0)  and  xa  €  R"  is  consistent  with  the  DAE  (3.17-)  (resp. 
(3.17+))  at  to,  the  initial  value  problem 

(3.18)  Ax  +  Bx  =  b,  *(to)  =  lo, 

has  a  unique  solution  x  e  Ci^p(R*).  Furthermore,  if  is  an  open  interval  about  to  and 
b*  €  Cj]),p(R*)  is  a  sequence  such  that  6*^  =  ,  for  all  C  >  1,  as  distributions  in  I,„  and 

lim  =  b  in  (2?'(R))",  then  xq  is  consistent  with  all  the  DAE’s  obtained  by  replacing  b 
by  6'  in  (3.17-)  (resp.  (3.17+),  and  denoting  by  x'  €  Cj^p(R*)  the  unique  solution  of  the 
initial  value  problem 

(3.19)  Ax'  -I-  Bx'  =  6',  x'(to)  =  xo, 
we  have 

(3.20)  limx'  =  x  in(Z>'(R))". 


IS 


(ii)  If  <0  =  0  wid  lo  G  is  consistent  with  the  DAE  (3.17^)  (resp.  (3.17^.))  at  to  =0, 
the  initial  value  problem 

(3.21±)  Ai  +  Di  =  5,  jr_(0)  =  lo  (resp.  i+(0)  =  i-fl), 

has  a  unique  solution  i  €  Furthermore,  if  Iq  is  an  open  interval  about  0  and 

6*  €  Cimp(R*)  is  *  sequence  such  that  5^^  =  6|,^,  for  all  f  >  1,  as  distributions  in  /o  and 

lim  b‘  =  b  in  (2)'(R))",  then  xo  is  consistent  with  all  the  DAE’s  obtained  by  replacing  b 
t— *00 

by  6^  in  (317-)  (resp.  (3.17+)]  and  denoting  by  x^  e  Cvmp(R‘)  ^fte  unique  solution  of  the 
initial  value  problem 

(3.22)  Ai'+Bi'  =  6*,xl(0)  =  io  (resp.  i'(0)  =  r„), 
ive  have 

(3.23)  Umx'  =  i  in  (D'(R))". 

Proo/;  In  this  proof,  we  use  the  notation  of  Section  2  without  further  mention.  From  the 
reduction  procedure  we  know  that  every  solution  x  €  (D'(R))"  of  the  DAE  (3.16)  has  the 
form  X  =  r^-ii,,  +  Ub-i  where  r„_i  and'ti„_i  are  given  by  (2.6)  and  (2.7),  respectively, 
and  i„  solves  the  ODE 

(3.24)  x„  +  BpX„  =  A„  ^by. 

The  key  point  here  is  the  simple  fact  that  the  distributions  Uj  and  b^  of  Section  2 
belong  to  Cj'iJ,p‘(R*)  and  have  impulse  order  at  most  k  +  j.  Indeed,  recall  that  bo  =  b 
and  r_i  =  n,  whence,  because  of  up  =  B^{AA^  +  BB^)  ~*b,  we  have  uq  £  and 

iord(uo)  <  iord(6)  =  k  by  Remark  3.2.  Therefore,  up  6  CjInp(R*)  has  impulse  order  at 
most  1:  +  1  which,  in  turn,  implies  that  bi  =  D(b  —  Buo  —  Auq)  6  CimpC®*)  impulse 
order  at  mostfc  +  1.  Obviously,  the  statement  about  the  sequences  up ,.. .  ,u„_i,6o,...  ,by 
now  follows  inductively  by  the  same  argument. 
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Since  6  has  impulse  order  at  most  t  +  —  1,  the  same  is  true  of 

Therefore,  by  Lemma  3.1,  the  solutions  of  (3.24)  are  in  ®nd  have  impulse  order 

at  most  it  +  i/  —  1.  This  implies  that  the  solutions  i  =  r„_)  +  tv  of  (3. 16)  are  in  C,^p(R*). 

Morover,  lord  (r^_ii„)  <  k  +  u  —  1  since  Fi,-!  is  C*,  and  iord  (tv)  5  L  +  i/  —  1,  because 
the  Cj's  are  C*'  and  iord  (u^)  <  k  +  j  for  0  <  y  <  t/  —  1.  Thus  t  has  impulse  order  at 
most  fc  +  r”  —  1 . 

If  now  to  <  0  (resp.  fo  >  0)  and  xa  6  R"  is  consistent  with  the  DAE  (3.17- )  (resp. 
(3.17+),  there  is  a  unique  io„  €  R'"'"'  such  that  lo  =  ry_i(<o)i./o  +  (note 

that  u,/-i(ta)  makes  sense  since  to  ^  0).  Hence,  the  solution  i  of  (3.18)  is  obtained  as 
X  =  +u„_i  where,  in  line  with  Theorem  3.1,  the  unique  solution 

of 

(3.25)  x^  +  A~' B^x„  =  b„,  x,(to)  =  i„o, 

and  no  other  initial  values  can  be  substituted  for  i„o  because  r,_i(to)  is  one-to-one. 

For  the  “furthermore”  part  in  (i)  of  the  theorem,  observe  first  that  consistency  of  a 
value  lo  €  R"  with  a  (linear)  DAE  at  a  point  to  depends  only  upon  the  coefficients  and 
the  right-hand  side  of  the  DAE  in  an  zubitrarily  small  neighborhood  of  to.  As  a  result, 
the  hypothesis  i^|,_  =  ensures  that  x'o  remains  consistent  with  the  DAE  obtained  by 
replacing  6  by  6*  in  (3.17_)  (resp,  (3.17+). 

For  fixed  f  >  1,  denote  by  uj,0  <  7  S  ~  1  tuid  bj,  0  <  j  <  i^,  the  sequences 
corresponding  to  Uj,  6,  in  the  procedure  of  Section  2  after  replacing  b  by  6*,  and  let  iij_, 
be  defined  by  (2.7)  with  uq,...  ,u„_i  replaced  by  uj,...  respectively.  With 

as  in  (2.6),  we  find  that  the  solution  i*  of  (3.19)  has  the  form 

(3.26)  =  + 

where  ij  €  solves  the  initial  value  problem 

(3.27)  ii  +  A:'B^x^  =  A:'bl  xiM  =  lo- 
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FVom  the  hypothesis  bf  =  ki,  ,  it  follows  at  once  that  ui  =  u„  and  k*  =  k,, 

I'-o  ^I'.o  •"'■o 

for  all  the  indices  f,  j  of  interest. In  particular,  k^  =  k„  and  hence 

^‘0 

(3.28)  A;'bi  =A:‘6  ,  W>1- 

''to  "to 

Next,  the  hypothesis  lint  k'  =  k  in  (jyfR))"  and  continuity  of  the  multiplication  of  distri- 
*00 

butions  by  C°“  matrix-valued  functions  yield  lim  tij  =  u,  in  (P'(R))rj_i,  0  <  j  <  i/  —  1 
and  ^lim  kj  =  bj  in  (I>'(R))'''-*,  0  <  j  <  u.  In  particular, 

(3.29)  lim  A;’bt  =  a;' k„  in  (P'fR))’'--' 

*00 

and 


(3.30)  limui.,  =u,_,  in  (P'{R))". 

t-*QC 

Since  i„  and  x*  solve  the  initial  value  problems  (3.25)  and  (3.27),  respectively,  it  follows 
from  (3.28)  and  (3.29)  and  Theorem  3.1  (iii)  that  lim  ij  =  Xy  in  (I>'(R ))'’’•-•.  Together 

/^OO 

with  (3.26)  and  (3.30),  this  implies  that  lim  i*  =  -f  v„_i  =  i  in  (I>’(R))".  This 

i-^OO 

completes  the  proof  of  part  (i)  of  the  theorem. 

Finally,  for  the  proof  of  (ii),  if  to  =  0  and  Xq  €  R"  is  consistent  with  the  DAE  (3.17- ) 
(resp.  (3.17+)  at  to  =  0,  then  the  solution  x  of  (3.21^)  is  obtained  in  the  form  i  = 
Ty-iXy  +  Ui,-i  where,  by  Theorem  3.1,  x„  is  the  unique  solution  of 

Xy  +  AZ'ByXy  =  by,  X (0)  =  X ,-0  (  resp .  X ,+ ( 0 )  =  r„o), 

and  i„o  €  R’’""'  is  (by  injectivity  of  r»-i(0))  the  unique  solution  of  the  equation  lo  = 
r„_i(0)x^o-f  u_(0)  (resp.  xg  =  r,,_i(0)i„o+*’+(0)).  The  proof  of  the  remaining  statement 
is  identical  to  the  proof  of  the  “furthermore”  part  in  (i)  of  the  theorem.  □ 

Remark  3.5:  Let  “(  ]"  stand  for  “jump  at  0”.  Then,  since  every  solution  i  of  the 
DAE  (3.16)  has  the  form  i  =  Vy-iXy  u„_i  with  x„  solving  the  ODE  (3.24),  we  have 
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[t]  =  r„-i(0)(i„]  +  and  limp  =  (r,_ii»)imp  +  IV  imp-  On  the  other  hand,  r^_i 

and  v„-x  are  obtained  through  an  explicit  procedure.  As  a  result,  [i]  and  ijmp  can  be 
calculated  if  [i„]  and  i„  imp  are  known  (using  (3.3)  -  (3.4)  for  the  term  (r„_ii^)in,p).  But 
from  Remark  3.4,  [i»]  and  imp  can  be  evaluated  from  b„  imp.  and  b„  imp  is  calculable 
since  6„  is  known  explicitly.  Thus,  both  [i|  and  limp  are  calculable,  at  least  in  principle.  □ 

4.  Inconsistent  initial  values. 

Let  A,B  €  C*'(R;£(R"))  and  6+  €  ^“([0,00); R")  be  given.  As  noted  in  the  Intro¬ 
duction,  the  problem  of  solving  the  initial  value  problem 


(4.1) 

A(t)x  +  B(t)i  =  6+(t), 

(4.2) 

i(0)  =  lo. 

for  arbitrary  lo  6  R"  that  is  not  necessarily  consistent  with  the  DAE  (4.1)  at  to  =  0,  often 
arises  when  a  known  function  on  (— oo,0|  verifies  i_(0)  =  Iq  and  is  to  be  extended 
into  a  solution  of  the  DAE  in  (4.1).  A  general  approach  is  suggested  by  the  following 
observation: 

Lemma  4.1.  Let  x-  €  C'“((-oo,0l;R'')  be  given,  and  suppose  that  xo  =  x-(0)  is 
consistent  with  the  DAE  (4.1)  at  to  ^  0.  Assume  further  that  the  pair  (A,B)  has  index 
1/  >  0  in  R  and  that,  in  the  notation  of  Section  2,  A„(t)  is  invertible  for  every  t  €  R.  Let 
1+  €  C^dO,  00);  R")  be  the  unique  solution  of  (4.1)  and  (4.2).  Set 

f  A(l)i-(f)  +  B(t)i_(t)  ift<0, 

(4.3)  6(t)  =  { 

U+(f)  ift>0, 

so  that  b  e  CjJ^p(R*)  (and  6imp  =  Oj.  Then,  the  function 

(  i_(t)  ift  <  0, 

(4.4)  x(t)  =  { 

U+(<)  ift>0. 
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verifies  I  €C"„p(R*)n(C«(R))"  and  is  the  unique  solution  o{  both  initial  value  problems 

(4.5)  A(  +  B(  =  b  inR.  ^-(0)  =  lo. 

(4.6)  Ai  +  B(  =  b  inR.  {+(0)  =  lo. 

Proof:  It  is  obvious  that  i  6  C|'J.,p(R*).  In  particular,  the  derivative  z  6  (^^(R))"  is  the 
function  given  by  i_(t)  for  t  <  0  and  by  x+(t)  for  t  >  0,  whence  Ax  +  Bx  =  fc  in  R  in 
the  sense  of  distributions.  By  definition  of  x_  and  ft,  Xq  =  x_(0)  is  consistent  with  the 
DAE  (3.17_)  at  to  =  0,  and  by  hypothesis  lo  is  also  consistent  with  the  DAE  (3.17+)  at 
to  =  0.  It  then  follows  from  Theorem  3.2  that  (4.5)  and  (4.6)  each  have  a  unique  solution 
in  C^p(R*).  In  both  cases,  this  solution  is  x  since  i_(0)  =  x+(0)  =  lo  by  continuity  of  x 
at  0.  □ 

Lemma  4.1  suggests  that  we  should  solve  (4.1)  for  inconsistent  lo  by  making  use  of  the 
extension  ft  of  ft+  in  (4.3).  This  approach  is  taken  in  the  following  result: 

Theorem  4.1.  Let  i_  6  C'“((-oo,0l;R")  and  ft+  €  C“(l0,oo);R")  be  given.  Suppose 
that  the  pair  (A,B)  has  index  t/  >  0  and  that,  in  the  notation  of  Section  2,  A^(t)  is 
invertible  for  every  t  €  R.  Then,  for  ft  €  CjJi,p(R*)  defined  by  (4.3),  there  exists  a  unique 
distribution  x  €  (D'(R))"  ivhich  solves 

(4.7)  Ax  +  Bx=:b  inR,  =  i-. 

Morever,  we  have 

(ij  X  €  CiJ^p(R*)  and  i  has  impulse  order  at  most  u  —  2. 

(ii)  x+  =  i||„  solves  the  DAE 

(4.8)  A(t)i+ +  fl(t)i+  =  ft+(t)  in(0,oo). 

(Hi)  If  xo  =  i-(0)  is  consistent  with  the  DAE  (4.8)  at  to  =  0,  then  x  6  (C°(R))”  and 
1+  is  the  classical  solution  of  the  initial  value  problem 

(4.9)  A(t)i+ +  B(t)i+ =  ft+(t)  in(0,oo),  i  +  (0)  =  xo. 
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(iv)  Irrespective  of  the  consistency  of  xt  =  i-(0)  with  the  DAE  (4.8)  af  (o  =  0,  the 
distribution  (  =  +  limp.  extended  by  0  in  (— oo,0).  is  the  unique  .solution  in 

C,mp(R*)  of  the  initial  value  problem 

(4.10)  Ai+  B(  =  b+  +  A(0)ioS  in  R,  {-(0)  =  0, 

where  6+  is  extended  by  0  in  (— oo,0).  In  particular,  iord  (^)  <  max(-l,  i/  —  2)  and  hence 
iord  (8)  <1^  —  2  for  >  1  (in  contrast  to  iord  (f )  <  u  —  1  obtained  by  a  direct  application 
of  Theorem  3.2  to  (4.10)). 

Proof:  Set  xo  =  x_(0)  and  6_(t)  =  j4(t)i_(t)  +  B(t)i-(<)  for  /  <  0,  so  that  by  definition 
x_  solves  the  DAE  (3.17_ )  and  xo  is  consistent  with  (3.17_)  at  to  =  0.  Thus,  by  Theorem 
3.2  there  is  a  unique  solution  y  €  <^imp(R*)  of  the  initial  value  problem 

(4.11)  Ay  +  By  =  b  in  R,  t/_(0)  =  xp. 

As  a  result,  y_  is  another  solution  of  (3.17_)  which,  just  as  x_,  verifies  !/-(0)  =  Xp.  This 
implies  that  y_  =  x_  and  hence  that  x  =  y  solves  (4.7).  Conversely,  if  y  €  (D'(R))"  solves 
(4.7),  then  y  €  Ci^p(R*)  by  Theorem  3.2,  and  the  equality  y|,.,„  =  x_  as  distributions 

in  (—00,0)  implies  at  once  that  y_  =  x_.  Thus,  by  continuity,  y-(0)  =  x_(0)  =  xp,  and  y 
solves  (4.11).  Uniqueness  of  the  solution  of  (4.7)  then  follows  from  the  unique  solvability 
of  (4.11). 

We  now  pass  to  the  proof  of  the  statements  (i)  -  (iv).  Part  (i)  follows  from  Theorem 
3.2  and  the  fact  that  the  right  side  6  in  (4.7)  and  (4.11)  (which,  as  was  just  seen,  have  the 
same  solution)  is  given  by  (4.3),  and  hence  has  impulse  order  it  <  —  1.  Property  (ii)  is  a 
trivial  consequence  of  (i)  and  the  fact  that  x  solves  (4.7).  For  the  proof  of  (iii)  note  that  if 
Xo  =  i-(0)  is  also  consistent  with  the  DAE  (4.8)  at  to  =  0,  then  Theorem  3.1  ensures  that 
the  unique  solution  y  =  x  of  (4.11)  (and  hence  also  of  (4.7))  is  given  by  (4.4)  and  solves 
(4.6).  This  shows  that  i  €  (C'’(R))"  and  that  x+  solves  (4.9). 

For  the  proof  of  (iv),  set  {  =  i+  +  Xj^p  €  CjJ^pfR*),  so  that  {  =  x  -  x_  with  x_ 
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extended  by  0  in  (0,oo).  Since  lord  (i)  <  i'  — 2  and  iord(i_)  <  —1,  we  find  that  iord(0  < 
max(  — 1,1/  —  2).  Moreover,  viewing  i_  as  a  function  of  <  €  R,  we  Imve  /li_  +  Bx.  = 
b-  —  AxoS,  where 

(  /l(<)i_(t)  +  B{t)x.(t)  if  <  <  0, 

I  0  if  <>  0. 

From  the  above  discussion  and  the  definition  of  b  in  (4.3),  it  follows  at  once  that  ,4f  +  = 

Ax  +  Bx  —  j4i_  —  Bx-  =  6+  -t-  ,4io^  Moreover,  we  have  Axai  =  .4(0)j'o^  ajxl  (-(0)  = 
i_(0)  —  x_(0)  =  0  whence  i  solves  (4.10)  and  thus  coincides  with  the  unique  solution  of 
that  problem.  Note  here  that  the  consistency  of  0  €  R"  with  the  DAE  A(t)j_  +  B(t)(_  = 
(6+)_(t)  =  0  in  (-00,0)  follows  from  the  fact  that  {_(<)  =  0  is  a  solution.  O 

Theorem  4.1  justifies  the  choice  of  the  solution  (,  of  (4.10)  (or  of  its  positive  part  f.^) 
to  represent  the  solution  i  of  (4.1)  when  xq  is  not  consistent  with  (4.1)  at  to  =  0.  Further 
justification  will  be  provided  by  Theorem  4.2  below.  For  the  time  being  observe  that  the 
characterization  (4.10)  of  {  =  x+  +  Xin,p  shows  that  the  extension  x  of  i.  as  a  solution  of 
(4.7)  depends  only  upon  lo  =  i-(0)  and  6+  and  hence  is  independent  of  i-(f)  for  t  <  0. 
For  the  case  when  A  and  B  are  constant,  the  characterization  (4. 10)  is  exactly  that  of 
[VLI<81],[Co82],  [G93],  but  the  equivalent  characterization  (4.7)  is  not  explicitly  noticed 
in  these  papers. 

It  is  noteworthy  that  for  index  1  problems,  f  solving  (4.10)  has  impulse  order  at  most 
—  1,  and  hence  (imp  =  0.  In  other  words,  (  is  a  function  with  a  possible  discontinuity 
at  the  origin.  A  simple  formula  can  be  given  for  the  jump  (+(0)  of  (  (recall  (_  =  0) 
without  using  the  more  cumbersome  general  procedure  outlined  in  Remark  3.5.  Indeed, 
the  relation  A(  +  B(  =  6+  +  A(0)xo4  reads  i4(+(0)5  +  A(d^/dt)  +  B(  =  i+  +  j4(0)xoi, 
where  d(/dt  denotes  the  usual  derivative  of  (  at  points  of  R*.  Clearly,  this  requires  that 
A(0)(+(0)  =  A(0)io  and  that  A(t)(<I(/dt)(t)  +  B(t)((t)  =  i+(t),  Vt  €  R*.  In  particular, 
for  <  >  0,  we  must  have  Qo{t)B(t)^^{t)  =  QQ(t)b^.(t)  where  Qo  €  C°°(R;£(R"))  is  as 
in  Section  2;  that  is,  <3o(0  projects  onto  a  complement  of  rge  zl(t),  (or  of  ext  rge  A{t)  in 
the  analytic  case).  By  continuity,  we  obtain  Qo(0)B(0)(+(0)  =  Qo(0)fc+(0).  Thus,  (+(0) 


22 


solves  the  system 


(4.12)  (.4(0)  +  <?o(0)5(0))4+(0)  =  A(0)io  +  <?o(0)5+(0). 

Conversely,  if  f+(0)  solves  (4.12)  then  i4(0){+(0)  =  >l(0)io  anJ  because  of  rge  .4(0)  n 
rge  Qo(0)  =  {0}  we  have  yo(0)B(0)f+(0)  =  ^o(0)6+(0).  It  turns  out  (see  [RR93a))  that 
invertibility  of  .4(f)  +  Qo(<)B(0  for  every  f  €  R  is  implied  by  the  index  1  assumption,  and 
hence  that  4+(0)  is  given  by 

?+(0)  =  (>1(0)  +  Qo(0)B(0)]-'(yl(0)ro  4^  yo(0)i.+  (0)). 


Note  that  {+(0)  =  lo  if  and  only  if  Qa(0)B.|.(0)xo  =  Co(0)6+(0)  (see  (4  12)),  a  condition 
that  is  easily  seen  to  be  equivalent  to  the  consistency  of  xq  with  the  DAE  (317+ )  at  to  =  0. 

In  Theorem  4.1,  the  function  (){f)  may  be  approximated,  in  the  sense  of  (I>'(R))",  by 
sequences  of  smooth  functions  b*  €  C‘®(R;R").  In  practice,  considering  such  a  sequence 
amounts  to  viewing  the  transition  from  i_  to  x+  as  the  limiting  case  of  a  perhaps  physically 
more  realistic,  situation  where  a  rapid  but  not  discontinuous  modification  of  the  input 
occurs  in  the  vicinity  of  f  =  0.  In  this  setting,  it  is  perfectly  reasonable  to  assume  that 
=  6_,  for  all  f  >  1,  in  some  interval  (-oo,  -a]  for  some  a  >  0  independent  of  (.  On  the 
other  hand,  the  function  x_,  has  a  unique  extension  as  a  solution  6  C“’(R;  R") 

of  the  DAE 

Ai^  +  Bx*  =  bt  in  R. 


In  fact,  I*  can  be  obtained  as  the  solution  of  the  initial  value  problem 


(4.13)  Ai*  +  Bx*  =  6/  in  R,  x*(to)  =  J’-(to). 

where  to  <  —a  is  arbitrarily  chosen.  Evidently,  it  would  be  desirable  that  the  sequence  i* 
tends  to  the  solution  x  of  (4.7)  in  some  sense.  That  this  is  indeed  true,  and  more  specifically 
that  lim  x*  =  i  in  (©'(R))"  follows  at  once  from  Theorem  3.2  and  the  hypotheses  b‘  =  b 

t~-*oo 

in  (— oo,  —a],  for  all  f  >  1,  and  lim  b*  =  bin  (I^fR))"  (just  choose  to  <  —a  in  (4.13)).  We 

t^OO 

record  this  result  in  the  following  form: 
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Theorem  4.2.  Let  i_  €  C^K-oo,  0),R)  an'i  6+  €  r’“°((0,oc );  R")  be  given  Suppose 
that  the  pair  (A,B)  has  index  v  >  0  and  that,  in  the  notation  of  Section  2,  A^it)  ‘s 
invertible  for  <  6  R.  Let  k  €  defined  by  (4  3)  and  let  b'  €  C’°(R;R"),  ^  >  1. 

be  a  sequence  such  that  b*  =  6_  in  ( -oo.  -a)  for  some  a  >  0  independent  oft  and  such  that 
lim  b^  =  b  in  (Z)'(R))”.  Denote  by  z*  €  C°^(R;R")  the  unique  extension  of  r_  as 

t—ac  '-oo-.! 

a  soJution  of  the  DAE 

Ai'  +  Bx'  =  b'  in  R, 

and  let  x  6  Cj'J„p(R*)  be  the  solution  of  (4.7).  Then,  we  have 

limi'  =  i  in  (D"R))". 


5.  Some  generalizations. 

Let  C  R  be  an  open  interval  and  let  S  =  (a,  ),gz  be  r;  nondecreasing  setjuenre  of 
points  of  RU  {±oo)  with  o,  <  a,+  i  if  either  a,  or  a,+  i  is  real  and  lim  a,  i  J .  Denote  by 

I  — ±» 

C,a,p(J  \  S)  the  subspace  o{V{J)  of  the  distributions  of  the  form  i  =  x  +  Xjmp  w-here  x  is 
a  function  such  that  ^  ^  n  ),  Vi  6  Z,  and  xjmp  is  a  distribution 

with  support  contained  in  5n^.  Equivalently,  if  S^,  is  the  Dirac  delta  distribution  at  a,, 
then  ijmp  'S  *  finite  or  infinite  linear  combination  of  derivatives  of  6,.  with  a,  €  J  With 
this  definition  of  Cin,p(,7  \  S),  we  have  C,mp(J)  =  if  J'  0  5  =  0. 

With  the  definition  Ci1„p(.7^\>5)  =  (Cnip(J/^\5))",  it  should  be  evident  how  to  formulate 
Theorems  3.1  and  3.2  for  the  case  h  €  Cj'J,p(  \5).  Of  course,  the  elements  of  C",„p(J  \  S) 
have  an  impulse  order  at  each  point  6  ,7  n  5.  It  may  only  be  useful  to  note  that  for 
a  given  arbitrary  sequence  €  R",  a  primitive  of  “ 

ffa, )  +  53*0  where  /fa,(<)  =  H(t  -  Oj)  if  a,  e  R.  H-ac{t)  =  1,  Hao{t)  =  0,  and  not 

II*-oc  hoi  Wa,  which  would  not  make  sense  when  53,3*-<x.  hoi  does  not  converge. 

For  0  6  it  should  be  equally  obvious  how  solutions  of  the  initial  value  problem 
Ax  +  Bx  =  b+  in  J+ =  J  ri{0,oo),  i(0)  =  lo. 
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can  be  defined  wlien  0  {  S,  6+  6  C|^|,U+\5),  and  To  €  R"  is  not  consistent  with  the  DAE 
Ai  +  Bx  =  6+  in  at  t  =0  For  problems  with  index  u  >  2  the  solutions  may  exliibit 
a  nonzero  impulse  at  to  =  0-  The  case  considered  in  Section  4  coiresponds  to  =  R, 
5  =  {oo). 

6.  Numerical  examples. 

For  index-one  problems  the  computation  of  the  jump  (4.12)  caused  by  inconsistent  input 
can  be  easily  incorporated  into  a  numerical  procedure  for  solving  the  initial  value  problems 
(4  1/2).  VVe  consider  here  a  recently  developed  solution  process,  (RR93b],  which  is  based 
on  the  reductif>n  procedure  of  (RR93a]  summarized  in  Section  2  above. 

Suppose  that  the  DAE  (4.1)  has  index  1.  For  agiven  step  h  >  0  set  1,  =  th.  t  =  0. 1 . 

and  consider  the  explicit  Euler  approximation 

(6.1)  A(f,)i(j,+i  -  Xi)  +  =  5+(t,). 

In  [RR93b|,  (Theorem  3),  it  was  shown  that  any  solution  io,ii,. . .  ,x,n  F  R"  verifies  for 
I  =  0, 1, . . .  ,  m  -  1  the  equations  Q{t,)B{t,)x,  =  Q(f,)fe+(f,'  and 

(6.2)  [A(t,)  -I-  Q(<.+,)5(<,+,))i,+,  =  (A(/,)  -  /iB(t,)]x,  +  /i6+(t,)  +  Q(t,+  ,  )6,(f,+  , ). 

Conversely,  for  sufficiently  small  h  and  any  given  lo  6  R"  such  that  Q(0)B(0)X|)  = 
Q(0)6+(0),  the  solution  Xo.ii,...  ,Xm  of  (6.2)  is  unique  and  solves  also  (6.1).  Smallness 
of  h  ensures  that  the  operator  .4(<,)  +  Q(ti+i  )  is  invertible,  given  that  invertibility 

of  A[t)  +  Q(t)B(t)  for  all  t  is  equivalent  to  the  index  1  assumption. 

The  difference  scheme  (6.2)  has  been  used  as  the  base  method  in  an  explicit  extrapola¬ 
tion  integrator,  LTVIXE,  for  general  index-one  problems  (4.1)  -  (4.2). 

Now  note  that  for  to  =  0,  /i  =  0,  and  with  i,+i  replaced  by  f+(0)  the  difference  equation 
(6  2)  is  identical  with  (4.12).  Thus,  the  results  of  Section  4  ensure  that  for  any  given  Xq 
we  only  need  to  apply  (6.2)  with  fi  =  0  to  obtain  the  consistent  starting  point  from  which 
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the  solution  process  can  then  be  started.  This  represents  only  a  minor  modification  to 
the  mentioned  code  LTVIXE.  The  resulting  code  accepts  any  given  initial  point  and  then 
computes  the  solution  starting  from  the  corresponding  solution  of  (4.12). 


As  an  example  consider  the  index-one  problem 


given  in  [CPS8]  which  has  the  general  solution 


(6.4) 


x(<)  =  (ate'  +  0e~',  ae'  +  tsinf,  sin/)'  6  R^. 


For  several  randomly  selected  points  tmd  starting  times,  Table  6.1  gives  the 

corresponding  consistent  starting  points  computed  by  LTVIXE.  It  is  readily  checked  that 
these  consistent  points  verify  (6  4)  for  suitable  constants  a  and  0. 


point-type 

/ 

Xt 

X2 

Xi 

given  pnt. 

0.0 

0.0 

1.0 

1.0 

0.0 

0.0 

1.0 

0.0 

given  pnt. 

1.0 

1.0. 

1.0 

1.0 

1.0 

1.0 

0.84147098 

0.84147098 

1.0 

-1.0 

4.0 

5.0 

1.0 

-1.0 

0.15852902 

0.84147098 

1  given  pnt. 

-5.0 

3.0 

2.0 

-5.0 

0.90929743 

-1.0 

1.0 

-1.0 

2.0 

-1.0 

1.0 

1.8414710 

-0.84147098 

Table  6.1:  Consistent  points  for  (6.3) 


For  the  index-two  case  a  code  LTV2XE  was  developed  which  incorporates  the  reduction 
discussed  in  Section  2  for  a  given  index-two  problem  (4.1)  and  then  applies  LTVTXE  to 
the  reduced  index-one  problem.  The  central  part  in  the  reduction  is  the  computation  of 
the  mappings  C  and  D.  This  can  be  implemented,  in  general,  by  using  a  singular  value 
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decomposition  (SVD)  to  obtain  a  basis  of  rge  A  and  the  projection  Q,  and  another  SVD 
for  generating  a  basis  of  ker  QB.  But,  it  turns  out  that  in  many  cases  there  are  much 
simpler  ways  of  generating  these  mappings.  Thus  LTV2XE  assumes  that  subroutines  are 
available  not  only  for  the  coefficients  A,  B,  b  and  their  derivatives,  but  also  for  C  2md  D. 
This  allows  us  to  bypass  eeisily  the  costly  general  method  for  calculating  these  matrices 
whenever  a  simpler  approach  is  feasible. 

In  this  form  LTV2XE  will  work  as  long  as  the  coefficients  of  the  problem  are  smooth. 
When  the  right  side  of  the  original  equation  has  a  jump,  then,  in  general,  the  right  side  of 
the  reduced  equation  exhibits  not  only  a  jump  but  also  an  impulse.  Hence  the  earlier  given 
simple  jump  computation  (4.12)  for  index  one  problems  is  insufficient  for  the  index-two 
case. 


As  an  illustration  consider  the  simple  DAE 


(6.5) 


'1  0  0\ 


^0  0  1' 


0  1  O  i-f  O  1  0  i  = 


.0  0  0/ 


a  0  o> 


with  the  initial  condition  1“  =  (0,0,1)^.  When  r  is  a  smooth  function  with  t(0)  =  1. 
f(0)  =  0  then  the  unique  solution  is 


(6.6)  ii(<)  =  T(t),  i2(<)  =  <  -  1  +  exp(-f),  X3(t)  =  1  -  f(<). 


Suppose  now  that  T(t)  =  Hi{t)  where  Hi  is  the  Heaviside  function  with  the  step  at  t  =  1. 
Then  the  solution  has  the  same  form  as  (6.6)  but  with  xi(t)  =  Hi(t)  and  T}(t)  =  1  -6i(<). 
Thus  at  t  =  1  we  have  a  jump  of  size  1  in  the  first  component  and  an  impulse  of  size  -1 
in  the  third  component.  A  graph  of  this  solution  does  not  show  the  impulse.  But  if  we 
approximate  the  step  of  Hi  by  a  cubic  spline;  that  is,  if  we  consider  (6.5)  with 


(6.7) 


r(t)  = 


0, 


for  0  <  <  <  1  —  e. 


I  -t-  i  <z(<)(3  -  <t(/)^),  for  1  -  e  <  <  <  1  -I-  e,  <?(<)  =  ifi 
1 ,  for  <  >  1  +  c. 
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with  small  e  >  0,  then  (6.6)  shows  that  ij(l)  =  1  —  3/(4e).  In  other  words,  this  solution 
approximates  the  impulse. 

For  the  general  computation  of  the  jump  and  impulse  in  the  index-two  case,  suppose 
that  for  the  given  DAE  (4.1)  we  have  6  =  5  +  [6]/f|„  where  6  is  the  smooth  part  of  the 
function  and  [6]  a  jump  at  the  time  to.  Then,  we  have  [uol^io  —  5(*o)^(A(to)A(/o)^ + 
fl(to)fl(fo)^)~’l^l^<o  which  implies  that  tio  has  the  impulse  B^(fu)(.4(/o)A(fu)^  +  B{to) 

S(fo)^)“*M^io-  Accordingly,  the  right  side  5i  =  D(b  -  Aiio  -  Buo)  has  the  impulse 

(6.8)  =  ~DAD'^(AA^  +  BB^)-'(to)l(>)^.- 

Now  let 

(6.9)  Ai  ii  +  Bi  X)  =  6i.  6,  =  6]  +  (6i]B,„  + 

be  the  reduced  equation  and  consider  its  solution  in  the  form  ij  =  i,  +  [ijIBij  +  {i6|„. 
By  substituting  this  into  (6.5)  and  comparing  terms  we  obtain  the  conditions 

■^i(*o)6=0,  ■4i(to)(3:il  +  BI(<o)^I  =  ^1,  Qi(<o)-Si(fo)li)l  =  f?i(<o)[l'il 

which  can  be  combined  into  the  two  systems 

(6.10a)  (Ai(fo) -f  Qi(fo)Bi(fo)){j  =  Qi(fo)^i 

(6.106)  (A,(fo)  +  Q.(to)B.(to))[i,l  =  <?.(to)16il  +  ^.  -B.(<o)ei. 

Since  the  DAE  is  assumed  to  have  index  two,  the  matrix  Ai(<o)  +  Qi(fo)Bi(to)  is  nonsin¬ 
gular  rmd  hence  the  two  systems  (6.10a/b)  can  be  solved  sticcessively.  A  brief  calculation 
shows  that  (6.10b)  reduces  to  (4.12)  exactly  if  —  Bi(fo)fi  =  0. 

The  relations  (6.8),  (6.10a/b)  were  incorporated  into  LTV2XE  to  adlow  for  the  compu¬ 
tation  of  the  jump  and  impulse  at  any  point  Ig  where  the  right  side  6  of  the  original  DAE 
(4.1)  has  a  jump. 
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As  numerical  example  we  consider  the  following  index- two  problem 


(6.11a) 


ii  1]  -  i2  -  ^4  -  =  6 

ij  +  I’l  -  +  txj  -  Is  =  0 

Xj  -  tZi  -  Xj  -  tX4  =  0 
X4+(t  -  1)X2  -1-  I3  -  <X4  =  0 

t^Il  -f  (1  -  t)^X2  -Kt  -  2)13  =  t(<) 


where  r(f)  =  1  for  t  <  1  and  r(f)  =  -1  for  f  >  1.  For  the  consistent  starting  point 


(6.116)  X|  =  0.5,  X2  =  0.0,  X3  =  -0.5,  X4  =  0.0,  xj  =  0.0, 


Table  6.2  shows  all  steps  computed  by  LTV2XE  for  the  problem  (6.11a/b)  for  0  <  t  < 
1.2.  A  relative  tolerance  of  10"®  and  a  maximal  step  of  0.1  was  used. 

The  discontinuity  at  t  =  1  causes  a  recalculation  of  the  point  obtained  at  <  =  1  from 
which  the  solution  proceeded.  Clearly,  in  order  to  capture  the  discontinuity  exactly  at 
t  =  1,  this  value  has  to  be  included  in  tire  list  of  required  output-points  of  the  code.  This 
was  indicated  in  the  table  by  a  dividing  line.  At  any  jump  point  the  output  of  the  code 
includes  the  values  of  the  jump  and  the  impulse  of  the  solution.  In  this  case,  we  found  that 
at  <  =  1  the  solution  has  the  jump  (— 2,  — 2,0, 0,2)^/f]  and  the  impulse  (0,0,0,0,2)^6i. 
Of  course,  the  jump  is  also  clearly  seen  in  Table  6.2. 

In  analogy  with  the  simple  problem  (6.5)  we  approximate  the  step  by 

11,  for  0  <  f  <  1  —  e, 

ja(f)®  —  |a{t),  for  1  -  <  <  <  <  1  -t  f,  <T(t)  =  ^  , 

—  1,  for  t  >  1  -P  e. 
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t 

^3 

t4 

1 

0.0000 

0.0000 

0.0000 

0.5000(-l) 

0.4724 

-0.2763(-l) 

0.2500(-l) 

0.3984 

-C.5751 

0.7502(1) 

0.2941 

-0.6263 

0.1252 

0.3475 

0.1612 

-0.3317 

-0.6788 

0.1752 

1.184 

0.4475 

0.3369(-2) 

-0.4816 

0.7384 

-1.351 

-0.1384 

-0.6120 

0.2875 

-1.193 

-0.2494 

-0.7057 

mmm 

-0.7034 

-0.2862 

-0.7284 

-0.9317 

-0.3105 

0.7919 

-0.2890 

-0.7094 

-1.003 

0.5225 

0.1445 

0.8639 

-0.2570 

-0.6571 

1.060 

-0.2073 

-0.5875 

-1.105 

EH  iBE 

0.9969 

-0.1452 

-0.5060 

-1.141 

0.8545 

0.5046 

1.000 

-0.1421 

-0.5020 

-1.142 

0.8608 

0.5020 

-1.142 

0.8608 

HlB 

-2.496 

-1.146 

0.8644 

2.504 

'S 

-2.493 

-1.149 

0.8663 

2.505 

1.010 

-2.111 

-2.468 

-1.167 

0.8818 

2.511 

1.034 

-2.387 

-1.224 

0.9341 

2.513 

1.058 

-1.281 

0.9910 

2.490 

1.082 

-1.887 

-2.220 

-1.334 

■Eh 

2.444 

1.105 

-1.813 

-2.136 

-1.386 

2.379 

1.129 

-1.740 

-2.053 

-1.435 

BIh 

2.296 

1.151 

-1.668 

-1.971 

-1.481 

2.197 

1.174 

-1.597 

-1.889 

-1.523 

2.086 

1.196 

-1.528 

-1.810 

-1.563 

BsH 

1.963 

1.200 

-1.516 

-1.797 

-1.569 

■m 

1.941 

Table  6.2:  Solution  of  (6.11a/b) 


with  small  e  >  0,  then  we  ex^jc  c:  the  solution  to  approximate  the  impulse  (0, 0,0,0,2)^ii. 
Figure  6.1  show  the  fifth  comp  vients  in  the  case  of  e  =  0.05.  For  smaller  values  of  e  the 
system  becomes  too  stiff  to  capture  the  impulse. 
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FIGURE  64 

It  may  he  noted  that  we  did  not  succeed  to  compute  the  solutions  of  this  problem  either 
with  DASSL  (see  e.g.[BCP89))  or  RADAU5  (see  (HW91J).  Both  codes  failed  at  start-up. 
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theoretical  results  yield  an  efficient  numerical  procedure  for  the  calculation  of  the  jump 
and  impulse  of  a  solution  at  a  point  of  discontinuity.  Numerical  examples  are  given. 
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